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ABSTRACT 

The 4 Ms Chandra Deep Field-South (CDF-S) and other deep X-ray surveys have been highly 
effective at selecting active galactic nuclei (AGN). However, cosmologically distant low-luminosity 
AGN (LLAGN) have remained a challenge to identify due to significant contribution from the host 
galaxy. We identify long-term X-ray variability (~month-years, observed frame) in 20 of 92 CDF-S 
galaxies spanning redshifts z m 0.08 — 1.02 that do not meet other AGN selection criteria. We 
show that the observed variability cannot be explained by X-ray binary populations or ultraluminous 
X-ray sources, so the variability is most likely caused by accretion onto a supermassive black hole. 
The variable galaxies are not heavily obscured in general, with a stacked effective power-law photon 
index of r sta ck 1-93 ± 0.13, and are therefore likely LLAGN. The LLAGN tend to lie a factor 
of «6-80 below the extrapolated linear variability-luminosity relation measured for luminous AGN. 
This may be explained by their lower accretion rates. Variability-independent black-hole mass and 
accretion-rate estimates for variable galaxies show that they sample a significantly different black- 
hole mass-accretion rate space, with masses a factor of 2.4 lower and accretion rates a factor of 22.5 
lower than variable luminous AGN at the same redshift. We find that an empirical model based on 
a universal broken power-law PSD function, where the break frequency depends on SMBH mass and 
accretion rate, roughly reproduces the shape, but not the normalization, of the variability-luminosity 
trends measured for variable galaxies and more luminous AGN. 
Subject headings: galaxies: active — X-rays: galaxies 



1. INTRODUCTION 

Observations show that all nearby galaxies with a 
massive bulge c omponent host superma s sive black holes 
(SMB Hs) (e.g., iFerrarese fc Fordl l200l IGultekin et all 
2009). SMBHs accreting near the Eddington limit 
(L/LEdd ~ 0.1-1) are visible as luminous active galac- 
tic nuclei (AGN) that often outshine their host galaxies. 
Models of AGN lifetime, constrained by observed Ed- 
dington ratio distributions, suggest that SMBH growth 
is dominated b y this luminous phase, lasting ~ a few x 
10 8 years (e.g. .'(Marconi et alj|200l iShankar et al.ll2004t 
iHopkins fc Hernquistll2009f ). 

Observations including constraints o n the sizes of 
ionized "bubbles" around quasars (e.g., [jakobsc n et al.l 
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I2003t iGoncalves et al.l |2008 D and the lengths of relativis- 
ticjets and radio lobes fe.g.. lScheuerll995tlBlundell et al.l 
1999) suggest that the episodic lifetime of luminous ac- 
tivity is similar to the total lifetime, implying that a 
SMBH is triggered to the luminous AGN phase no more 
than a few times. SMBHs therefore spend significant 
amounts of time in quiescent or low-activity phases, 
which may contribute up to ~20% of overall SMBH 
growth (|Hopkins fc HernquistJ 12009ft . A SMBH accret- 
ing at lower rates (L/LEdd *C 0.1) will appear as a 
low-luminosity AGN (LLAGN). LLAGN share several 
properties with more lumi nous AGN, includin g similar 
X-ray spectral shapes (e.g.. lYounes et afll2011l) and sim- 
ilar radio-loud fractions and luminosit y-dependent, in- 
trinsic X-ray to optical flux ratios (e.g.,[Maoz 200jj). A 
more complete census of LLAGN is important for under- 
standing SMBH accretion history, but the relative signif- 
icance of the host galaxy in LLAGN makes a full census 
of accretion activity a challenge. 

Deep X-ray surveys have been effective at selecting a 
wide variety of AGN, including luminous, unobscured 
AGN as well as faint an d/or obscured AGN (e.g., see 
iBrandt fc Hasingcr 2005, for a review). X-ray selection 
criteria usually include cuts on X-ray luminosity and 
X-ray spectral shape. Multi-wavelength data further aid 
X-ray selection by allowing selection via, for example, ex- 
cess X-ray emissi on compared to what is ex pected from 
optical flux fe- g.-lHornschemeier et af1l2003[ ) or radio lu- 
minosity (e.g lAlexander et al.ll2005F iT 

The above methods have been successful in select- 
ing a wide variety of AGN, but nevertheless miss cer- 
tain populations, s uch as very heavily obscured AGN 
and LLAGN (e.g., iBauer et all 120041 : lAlexander et al.1 
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[200llLehmer et al.ll2008l ). While hea vily obscured AGN 
can often be selected in the IR (e.g.. iHouck et al1l2005l : 
iMartinez-Sansigre et all 120051: lAlexander et al.l 120081) . 
the spectral energy distributions (SEDs) of LLAGN are 
likely dominated by the host galaxy in other bands. Even 
in X-rays, X-ray binaries, ultraluminous X-ray sources, 
and hot gas will provide significant contributions to the 
overall power output. Simulated Chandra observations 
of nearby low-luminosity Seyfert nuclei artificially shifted 
to z ~ 0.3 show that LLAGN would exhibit X-ray lumi- 
nosities, spectral shapes, and X-ray-to-optical flux ratios 
consistent wit h those of normal or o ptically bright /X-ray 
faint galaxies ([Peterson et al.l [20061 ) . By relying on such 
criteria, deep X-ray surveys may be underestimating 
AGN fractions. 

Variability potentially provides a useful indicator of 
whether an extragalactic X-ray source, classified as a 
galaxy by other means, harbors an AGN. Variability is 
a defining characteristic of AGN an d has long been used 
as an AGN selection technique (e.g. . Ivan den Bergh et al.l 
1973). Numerous studies have used optical variability to 
select AGN from deep surveys such as the 1 Ms CDF-S, 
the Subaru/XMM-Newton Deep Fi eld, and the GOODS 
North and Sou t h Fields (e.g. . iTrevese et al.l 120081: 
Morokumal [ 200 8: Villfo rth et all 120101: iSaraiedini et all 



20111 respe ctively). Spectroscop ic observations of the 1 
Ms CDF-S (|Boutsia et all 120091) found that 17 of 27 op- 
tical variability-selected objects were broad-line AGN; 9 
(5) AGN would have been missed if selected by color 
(X-ray selection). 

Similarly, UV variability has been used successfully 
to identify LLAGN in galaxies with low-ionization nu- 
clear emission-line regions (LINERs). LINERs have been 
found i n the nuclei of a large fraction of nea rby galax- 
ies (e.g.. lHo et al.l[l997t iKauffman et al.ll2003l ). but these 
regions could be ionized by either massive star clusters 
or low accretion-rate AGN. HST imaging has found that 
~25% of LINE Rs are associated with compact (< few 
pc) UV sources (|Maoz et al.lll995l : iBarth et al.lll998D . A 
study of LINERs with compact nuclear UV sources found 
significant v ariability in 15 of 17, indicating the presence 
of an AGN (jMaoz et al.ll2005l ). 

Deep X-ray surveys are able to detect variabil- 
ity in moderate- luminosity /high-redshift AGN (e.g. . 
Almai ni et aT1l2000l; iPaolillo et aLll2004i: iPapadakis et al 



The 4 Ms Chandra Deep Field-South (|Xue et al 



20111 ). the deepest X-ray survey to date, allows a prelim- 



inary classification of AGN on the basis of several ob- 
served quantities (see Sj5]for details). This paper utilizes 
X-ray variability techniques to search for AGN missed 
by these criteria. With 4 Ms of exposure time spanning 
10.8 years for 466 good-quality sources (see $3|), variabil- 
ity can be detected in sources with time-averaged fluxes 

" x . We 
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as faint as -Fo.5-8kcV ~ 5 x 10 ergs cm" 
use a cosmology with H = 70.4 km s" 1 Mpc" % i 
0.272, and 0\ = 0.728 (e.g.. IKomatsu et al.ll2011j) . 



2. OVERVIEW OF THE 4 MS CDF-S CATALOG 

The details of th e 4 Ms CDF-S source catalog are avail- 
able in iXue et al.1 ([20111 ): we provide a brief summary 
here. The 4 Ms CDF-S, constructed from 54 Chan- 
dra observations over 10.8 years, covers an area of 464.5 
arcmin 2 and reaches highest sensitivities of i*b. 5-2keV ~ 
9.1 x 10~ 18 ergs cm -2 s^ 1 and F 2 _ 8 k c v ~ 5.5 x 



10 -17 ergs cm -2 s _1 , with multi-wavelength coverage in 
more than 40 bands from the radio to the UV. Source 
candidates are detected using a 10 ~ 5 false-positive prob- 
ability threshold in wavdetect ([Freeman et al.1 [2002] ) 
and are then pruned using a binomial no-source prob- 
ability (see Appendix A2 of Weisskopf et al. 2007) P 
< 0.004 to obtain a more conservative list of 740 main- 
catalog sources, all of which are consistent with being 
point sources. Source extraction and photometry were 
conducted with ACIS extract (AE; Broos et al. 2010). 
AE models Chandra's High Resolution Mirror Assem- 
bly point spread function (PSF) using the MARX ray- 
tracing simulator^ The PSF model is used to generate a 
polygonal extraction region for each source that approx- 
imates the ~90% encircled energy fraction (EEF) con- 
tour of a local PSF measured at 1.497 keV. For <6% of 
the candidates, the sources are crowded (i.e., the polyg- 
onal source regions overlap) and smaller extraction re- 
gions that are as large as possible without overlapping 
(40 - 75% EEF) are used. The background is calcu- 
lated from regions that subtract the contribution from 
the source of interest and its neighboring sources; the 
regions are typically a factor «16 larger than the source- 
extraction region. AE merges the individual observations 
to estimate aperture-corrected, background-subtracted 
counts and th e la (asymmm etric) upper and lower sta- 
tistical errors (jGehrelsill986l ). In this paper, we will use 
the standard X-ray photometric bands: 0.5-2 keV (soft), 
2-8 keV (hard), and 0.5-8 keV (full). 

Though most sources have a relatively small number 
of counts (median net counts ss 77), a rough estimate of 
source spectral shape can be made by relating the band 
ratio (i.e., the ratio of the count rates in the 2-8 keV 
and 0.5-2 keV bands) to an effective power-law photon 
index, T e ff (F v oc v~ a = v~ r+1 ). For low-count sources 
where r o ff cannot be determined reliably, r e g is set to 
1.4, the sta cked (co-added) s pectrum of all sources in 
the CDF-S ([Tozzi et alJl200lKXue et al.ll2Tji~l[ ). which is 
consistent with th e unresolved spectrum of th e cosmic 
X-ray background ([Hickox fc Marke vitch 2006|). 

Of 740 X-ray sources, 716 (96.8%) contain matches 
in at least one of seven optical/near- infrared/radio 
(ONI R) catalogs: (1) t he ES Q 2.2 m WFI i?-band cat- 
alog (IGiavalisco et all [200l . (2) the GOODS-S Hub- 
ble Space Telescope (H ST) version r2.0z z-band catalog 
(IGiavalisco et alj200l (3) the GEMS HST z-band cata- 
log (ICaldwefiet alll2008l) (4) the GOODS-S MUSIC cat- 
alog (IGrazian et al.lT2006h . (5) the MUSYC iVband cat- 
alog (ITavlor et al.ll2009fl. (6U he SIMPLE Spitzer/IRAC 
3.6um catalog (iDamen et al.ll201lD. an d (7) the VLA 1.4 
GHz radio catalog ([Miller et all |2008[ ). 

Of 716 X-ray sources with multi-wavelength iden- 
tifications, 419 (58.5%) have spectroscopic redshift 
measurements , coll e cted from ILe Fevre et all (120041) . 
ISzokolv et all (12001 . IZheng et all (120041) IMignoli et al l 
(120051) , iRavikumar et all (120071). iVanzella et all 
Bom. iPopesso et all (I2009D . iTreister et all T2 009) , 
iBalestra et al.l ([2010D . and ISilverman et all ([2010D . A 
total of 343 (81.9%) of the 419 spectroscopic redshift 
measurements are "secure", in that they are measured 
at > 95% confidence levels with multiple secure spec- 



12 MARX is available at |http://space. mit.edu/CXC/MARX/index.htmll 
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tral features. 668 (93.3%) sources have high-quality, 
accurate (|Az|/(l + z) » 6.5%) photometric-redshift 
measurements fro m at least one of t hree p h otometric- 
redshi ft catal ogs: iCardamone et all ()2010h . iLuo et ahl 
(|2OT0t) , and IRaffertv et all (|2011l ). The positions of 
primary ONIR counterparts were cross-matched with 
the photometric-redshift catalogs using a matching 
radius of 0.5", resulting in a false- match probability of 
<1%. Subsequent spectroscopic observations published 
in the Arizona CDF-S Environment Survey (ACES; 
Cooper et al. 2012) catalog show with a blind test that 
errors on the photometric redshifts are <1%. 

The 4 Ms CDF-S X-ray sources were classified as AGN 
by the following criteria: 

• High luminosities: Lo.5-8koV > 3 x 10 42 ergs s -1 , 
where the rest-frame luminosity has been corrected for 
Galactic and intrinsic absorption. 

• Hard spectra: A source with r c ff < 1.0 is identified 
as a heavily obscured AGN. 

• High X-ray-to-optical flux ratios: log(Fx- / Fr) > — 1, 

where F x = fb.s-skeV, -Fo.5-2kcV, or iVskeV and F R is 
the i?-band flux. 

• Excess X-ray emission compared to that expected 
from star format i on: L n.R-sksV > 3x(8.9 x 10 17 Lr) 
(|Alexander et al.l [2005), where L R is the 1.4 GHz 
monochromatic luminosity in W Hz -1 . 

• An indication of broad emission lines in the optical 
spectrum. 

Stars were classified by cross-matching X-ray sources 
(using the ONIR counterpart po sitions) with ( 1) th e 
spectroscopica l ly ide ntified stars inlSzokolv et al.l (12 004). 
iMignoli et aD (pOOl . and iSilverman et al.l (|2010h :72) 
the likely stars with stellarit y indices greater tha n 0.7 
in the GEMS HST catalog ((Caldwell et al.l 120081) : and 
(3) the likely stars with best-fit stella r templates in the 
MUS YC photometric-redshift catalog (jCardamone et al.l 

sou, 

using a matching radius of 0.5". 

X-ray sources not identified as an AGN or a star were 
classified as galaxies. 

Rest-frame 0.5-8 keV luminosities are calculated for 
all sources. For AGN, which make up the vast major- 
ity of CDF-S sources, the luminosity is corrected for 
Galactic and intrinsic absorption. AGN not detected 
in the full band have an upper limit on the X-ray lu- 
minosity based on the 3er Poisson error on the counts. 
For galaxies, the intrinsic-absorption correction and K- 
corrections may not be appropriate. Of 92 CDF-S galax- 
ies with "good quality" observations (see $3]) , 78 are not 
detected in the hard band and have poorly determined 
photon indices. All galaxies are detected in the soft- 
band. We perform a stacking analysis for galaxies with 
< 150 net counts (87 of 92) following the procedure de- 
scribed in ij5.ll The resulting average r sta ck = 1.90±0.08 
gives Lo.5-8kcv/^o.5-2kcV = 2.35. We calculate the 
rest-frame 0.5-8 keV luminosity for CDF-S galaxies as 
£o.5-8kcV = 47rd| x 2.35 x / .5-2kcv(l + ^) ratack " 2 - Since 
galaxy X-ray emission is typically unabsorbed, we do not 
apply any correction for intrinsic absorption. 



3. TESTING FOR X-RAY VARIABILITY 

We perform two quality cuts before conducting vari- 
ability analysis. First, we exclude the catalog sources 
with off- axis angles greater than 8' to ensure that sources 
will have sufficient coverage (>50 of 54 observations) 



throughout the various CDF-S pointings. To ensure ac- 
curate variability measurements, we also require that 
each source has at least 20 net counts in the 0.5-8.0 keV 
band (i.e., at least 5 counts on average per epoch, as de- 
fined below). These quality cuts result in a total sample 
of 466 CDF-S sources: 369 classified as AGN 92 classi- 
fied as galaxies, and 5 classified as stars in the lXue et al.l 
(|2011l) catalog. All of these AGN and galaxies have mea- 
sured spectroscopic or photometric redshifts. The 92 
sources classified as galaxies may nevertheless contain 
an AGN, as would be indicated by X-ray variability, and 
make up the sample investigated in this paper. 

We divide the CDF-S observations into four epochs, 
each containing ~1 Ms of integration: 2000 (943.1 ks), 
2007 (967.7 ks), 2010a (March-May ; 1015.5 k s ), an d 
2010b (May- July; 944.9 ks). As in IXue et al.l (|20rl . 
we merged observations within each epoch and, for a 
given source position from the CDF-S catalog, measured 
the source and background counts and flux over three 
observed-frame energy bands: 0.5-8 keV, 0.5-2 keV, 
and 2-8 keV. A source is considered variable if the 
variability observed between observations is greater than 
that expected from Poisson statistics, with a probability 
threshold of 5% that the variability is due to noise alone. 
(The choice of probability threshold is discussed further 
below.) To check whether a source can be considered 
variable, we calculate the quantity: 



N , 

E\ x i 

i=l 



(1) 



where N — 4 is the number of epochs, Xi is the pho- 
ton flux (background-subtracted counts with units of 
cm~ 2 s _1 ) in a given epoch, \i is the mean photon flux 
over all epochs, and a 2 is the error squared on the pho- 
ton flux for the i th epoch. The photon flux is calcu- 
lated by dividing the full-band (0.5-8.0 keV) net counts 
by the exposure time and the mean effec tive area across 
the source aperture. The iGehrel s (1986) approximation 
gives the error on the net counts, which is propagated 
to obtain the error on the photon flux. Since this error 
is significantly asymmetric for low-count (< 15 counts) 
sources, we average the upward and downward error bars 
for these objects to obtain the average error ct;. (The 
same method is applied in the Monte Carlo simulations 
below.) 

For large photon fluxes, the X 2 statistic follows a \ 2 
distribution, and any source with X 2 > 7.82 (for 3 de- 
grees of freedom) has a probability P X 2 < 0.05 (i.e., 95% 
confidence level) that the variability is due to random 
noise. However, at low count rates, the error on the 
photon flux is not Gaussian. Since errors in the low- 
count regime are larger than expected from a Gaussian 
distribution, the resulting X 2 statistic is smaller than ex- 
pect ed and does not fol low the x 2 distribution (see Fig. 
2 of lPaolillo et"aT1l200l . 

We therefore constructed a Monte Carlo simulation to 
determine the distribution the X 2 statistic should fol- 
low for each s o urce, similar to the procedure followed by 
iPaolillo et"~aT1 (|2004l ). We first scaled the total observed 
source and background counts for each source, obtained 
from the full 4 Ms observation, to the exposure time and 
effective area for a given epoch. This procedure gener- 



4 



ates the source and background counts expected in each 
epoch if the source and background were constant over 
time, and it accounts for fluctuations in the background 
that will affect low-count sources. To simulate the vari- 
ance expected from noise, Poisson distributions were de- 
fined using the expected source and background counts 
as the mean values. We then simulated 1,000 observa- 
tions of each source by repeatedly drawing the expected 
counts from the Poisson distributions. For each simu- 
lated observation, we calculated the photon flux for four 
epochs and calculated X 2 as defined above. Asymmetric 
error s on the s ource and background counts are obtained 
from IGehrelsl (fl98(l and are propagated to get the er- 
ror on the photon flux. The observed X 2 can then be 
compared to the simulated distribution to determine the 
probability Px 2 that the observed variability is due to 
Poisson noise. Spurious sources of variability are negli- 
gible, since effective exposure maps are calculated sep- 
arately for each observation, taking into account issues 
such as vignetting, CCD gaps, bad pixels, bad columns, 
and Chandra's spatial- and time-dependent quantum ef- 
ficiency degradation. 

While using a more conservative P cr it = 1% on our 
dataset would result in fewer false positives, it would also 
eliminate a similar number of truly varying sources. For 
example, in the sample of 92 galaxies, P cr it = 5% results 
in 20 variable sources (see S|4|), of which 4.6 are expected 
to be false positives. Reducing P cr ;t to 1% results in 13 
variable sources, of which 0.9 is expected to be a false 
positive. So while the more conservative critical value 
eliminates «4 false positives, it also eliminates w 3 — 4 
truly varying sources. The Px 2 values listed in Table 2 
can be used to screen the sources further as desired. 

4. GALAXIES WITH AGN-LIKE VARIABILITY 

Of the 369 CDF-S sources classified as AGN that meet 
both the off-axis angle (9 < 8') and count (total net 
counts > 20) requirements, 50.1% exhibit significant flux 
variability (Px 2 < 0.05) on ^month-year timescales. For 
the 178 AGN with more than 100 counts, 74.2% exhibit 
significant variability. The basic diagnostic plot in Fig- 
ure [1] shows the rise in the AGN variable fraction with 
total net counts. The plot demonstrates that, given suf- 
ficient counts to detect it, variability is a near- ubiquitous 
property of faint, X-ray selected AGN, even in the case 
of significant obscuration: «70% of the CD F-S sample 
consists of obscured AGN (|Xue et all 120111 ). and 47.5% 
(51.5%) of obscured (unobscured) AGN are significantly 
variable. The variable fractions are consistent wit h the 
results from the 1 Ms CDF-S (jPaolillo et all 120041 ). al- 
though obtained with different temporal sampling and 
down to much fainter fluxes. 

In sources classified as galaxies, the variable fraction 
is significant at low counts and equals that of AGN at 
higher counts (Figure [1]). Table 1 describes the columns 
of Table 2, which lists the attributes of the 20 variable 
and 72 non- variable galaxies. (Variability properties of 
CDF-S AGN will be covered in the forthcoming Pao- 
lillo et al, in preparation.) The luminosity-redshift dis- 
tribution of variable and non- variable galaxies is shown 
in Figure [5] almost all lie below z ~ 1. Spectroscopic 
redshifts are available for 18 of 20 variable galaxies and 
for 61 of 72 non- variable galaxies. Photometric redshifts 
are available for the remaining sources. Six example 



light curves (background-subtracted count rates in the 
observed-frame 0.5-8 keV band vs. MJD) representative 
of the sample as a whole are shown in Figure [3l 

The X-ray luminosity distributions of all CDF-S galax- 
ies and those exhibiting significant variability are shown 
in Figure |U A K-S test shows that the two samples 
are consistent with being drawn from the same parent 
population (Pks = 57%), and the variable fraction does 
not show a significant dependence on X-ray luminosity 
below i .5-8kov ~ 10 43 ergs s _1 . The possibility that 
AGN-relatcd variability may go undetected in galaxies is 
discussed in §4.11 

We briefly compare AGN selection based on variability 
to the following selection methods: (1) X-ray luminosity 
cuts, (2) the X-ray-to-optical flux ratio, and (3) excess 
X-ray emission compared to that expected from star for- 
mation, based on the radio luminosity. Figure [5] shows 
the fraction of variable sources vs. X-ray luminosity. Be- 
low io.5-8koV = 10 ergs s , a luminosity cut often 
used for AGN selection in X-ray surveys, the variable 
fraction remains significant at 20 — 30%. Of 64 CDF-S 
galaxies with L .5_8kcV < 10 42 ergs s _1 , 17 («27%) are 
variable. 

AGN selection via the X-ray-to-optical flux ratio is 
demonstrated by the i?-band magnitude vs. X-ray 
flux plane in Figure [5] (cf. Figure 16 in Xue et al. 
20 11). Sourc e s clas sified as AGN, galaxies, and stars 
in iXue et all (|2Qlll ) are marked as small red circles, 
larger black circles and blue stars, respectively. Variable 
sources are marked with filled symbols. The optically 
bright, X-ray faint region, typically thought to exclude 
AGN (OBXF; F X /F R < -2; Hornschemeier et al. 2003) 
contains 27 galaxies. Of these, 6 (22%) are variable, with 
X-ray luminosities spanning log Lo.5-8kcV ~ 39.7-41.4. 

AGN may also be selected based on a comparison be - 
tween X-ray and radio luminosities (|Xue et al.l 120 111 ) . 
The radio luminosity can be used to pre dict the X-ray lu- 
mino sity in star-forming galaxies fe.g.. lAlexander et all 
120051 and references therein), so sources with excess 
X-ray emission may be classified as AGN. Of the 17 
CDF-S galaxies with radio detections (none of which 
have excess X-ray emission), three (18%) are variable. 
These objects may have excess radio emission due to 
strong radio cores rather than star formation. 

4.1. Measuring Variability Strength 

Due to the generally limited photon statistics of the 
CDF-S galaxy sample, most variable sources must be 
strongly variable to be detected. Significantly vari- 
able galaxies show maximum-to-minimum flux ratios 
^max/min ~ 1-5-9.3 (median = 4.1) over the observed 
10.8-year time frame. The smallest max-to-min ra- 
tio (1.5) was measured for 033246.77-274212.7 (XID = 
616), a source with >500 counts. For most galaxies, to- 
tal net counts are too low (< 100) to detect variability 
below a factor of 2-3. 

To address whether the variable galaxy max-to-min ra- 
tios sample the average AGN population or only the 
highly variable "tip of the iceberg," we ran a Monte 
Carlo simulation. The procedure assumed that the en- 
tire galaxy population is significantly variable and simu- 
lated the variability expected ove r 10.8 years of observa - 
tion. Follo wing the procedure in iVaughan et all ((2003), 
we used the iTimmer fc Koenlgl (|1995l ) algorithm to sim- 
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ulate 5,000 light curves based on the mean and variance 
of the flux for each of the 92 CDF-S galaxies. For non- 
variable galaxies, the measured variance represents noise 
in the measured light curve, which gives an upper limit to 
the variability that could be present. The algorithm pro- 
duces a random, continuously sampled light curve from 
a given power spectral density (PSD) function, which we 
assumed to be described by a broken power-law, where 
the break frequency depends on SMBH mass and accre- 
tion rate (jMcHardv et all [200fih . Since the break fre- 
quency lies outside the ran ge of timescales sampled for 
most variable galaxies (see §6.3j) . we simplify the model 
to a power-law (P(f) oc J* 3 ) with index j3 — —1, as is 
typical for the low-fre quency (long-timesca le) PSDs for 
nearby Seyferts (e.g.. iVaughan et aD l2003h . The time 
units of the light curves are determined by the mini- 
mum and maximum timescales input into the simulation; 
the light curve durations are adjusted according to each 
source's redshift. 

We resampled/rebinned the light curve using the 
CDF-S observing pattern and add Poisson noise to the 
simulated light curve to account for measurement error. 
The full simulated light curves was made five times longer 
than the sampled region in order to produce variation on 
timescales much longer than those sampled by the data. 
This reproduces the effect that very long-timescale (low- 
frequency) variations have on variability measured over 
shorter timescales (i.e., "red noise leak"). 

We compare the simulated distribution of median max- 
to-min ratios, where the median is calculated over 5,000 
trials for each galaxy, to the observed distribution for 
variable galaxies in Figure [7] A K-S test shows that these 
two populations have a 0.2% chance of being drawn from 
the same parent population. Note that the simulated 
distribution (black histogram) illustrates an upper limit 
to the variability that could be present in the CDF-S, 
resulting in a lower limit on the detectable fraction of 
sources. The CDF-S detects at least «18% of sources 
with max-to-min ratios > 2 and at least «57% of sources 
with max-to-min ratios > 4. A significant fraction of 
non-variable galaxies may still host an AGN, but the 
variability may remain undetected due to low counts. 

For the galaxies exhibiting significant variability, 
we calculated the normalized excess variance (e.g., 
IVaughan et al.l l2003f) , which measures how strongly a 
source varies in excess of the measurement error. The 
excess variance is the integral of a source's PSD function 
over a given frequency range, which is defined by the 
light curve's duration (10.8 years, observed-frame) and 
minimum bin size (4.0 MsF^I observed- frame). 
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where a er r,i is the average of the asymmetric upward 
and downward measurement errors. Using the upward 
(downward) error would overestimate (underestimate) 
the errors. Zero excess variance (cr 2 xs = 0) would indi- 
cate that the observed count fluctuation is entirely con- 

13 While each CDF-S epoch totals ~1 Ms in integration time, 
the timescale sampled is significantly longer due to the spread in 
Chand ra observations, of which the shortest is 4.0 Ms for the 2007 
epoch IILuo et al.l(200D . 



sistent with noise rather than due to intrinsic source vari- 
ability; due to statistical fluctuations, the excess variance 
may also be negative in this case. The <7 2 xs values are 
listed in Table 2. Note that the variability amplitude is 
calculated for observed- frame energy bands. The vari- 
able galaxies cover a redshift range z = — 1, so the 
excess variance will be measured over 0.5 — 8 keV at z = 
to 1 — 16 keV at z = 1. Variability strength is known to 
depend on energy in some nearby Seyfert galaxies (e.g., 
Ark 120, MCG-6-30-15, and I Zw 1; Vaughan et al. 2004, 
Vaughan & Fabian 2004, and Gallo et al. 2007, respec- 
tively), while in others, variability remains nearly con- 
stant with energy (e.g., Ton S180 and NAB 0205+024; 
Vaughan et al. 2002 and Gallo et al. 2004, respectively). 
In the former cases, the change in variability strength is 
small, with a < 10% difference from 0.5 to 10 keV, so 
the bandpass effects at redshifts z — — 1 should remain 
small. 

We calculate the statistical error (i.e., measurement 
error) on the excess variance following Equation 11 of 
IVaughan et~aTl (|2003D . 
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The large errors on <r 2 xs (Table 2) are due to the small 
numbers of counts observed for most sources (e.g., 12 of 
20 variable galaxies have < 100 net counts); four vari- 
able galaxies with net counts < 50 have excess variance 
measurements completely dominated by statistical un- 
certainty [a„ xs < err(c 2 xs )]. Nevertheless, most variable 
galaxies have excess variance measured at the >lcr level. 

The excess variance contains additional sources of un- 
certainty aside from statistical error: (1) random scat- 
ter intrinsic to the stoch astic nature of AGN variabil- 
ity (jVaughan et al.1 120031 ) and (2) uncertainty and sys- 
tematic bias due to sparse sampling of the light curve. 
The sparse pattern of CDF-S observations will lead to 
large uncertainty in the mean flux measurement, and 
since the measured mean will be closer to the sampled 
data points rather than the true mean, th e variance mea- 
surements will tend to be underestimated (jAllevato et al.l 
lin preparation! ). We again employ a Monte Carlo simu- 
lation to model the uncertainty and systematic effects 
produced by random scatter and sparse sampling. 

We follow the same procedure described above to pro- 
duce 5,000 light curves for every variable source. The 
mean and variance are calculated after sampling each 
simulated light curve with the CDF-S observing pattern. 
The sampling bias can be corrected by rescaling the ob- 
served variance by a factor equal to the ratio between the 
"true" input variance (i.e., the observed variance used as 
input in the simulations) and the median output vari- 
ance (i.e., the biased variance affected by sparse sam- 
pling): / sca i c = cr I 2 ncas /median(o'2. m ). The median output 
variance is calculated over all 5,000 light curves. The 
amount of systematic bias depends on the frequency and 
regularity of the sampling. If the sampling is regular, 
the scaling factor will approach unity as the number of 
samples increases; however, the scaling factor will remain 
above unity even a t high sampling frequency if th e sam- 
pling is irregular (|Allevato et alJIin preparation! ) . The 
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slope of the PSD will also affect the sampling bias - 
a steeper PSD slope (i.e., (3 = — 2 instead of /3 = —1) 
will result in a larger bias for a given sampling pattern. 
Since the intrinsic PSD slopes are not known, we apply 
uniform corrections assuming j3 = — 1; however, source- 
to-source variations in PSD slope may be a significant 
source of scatter in variability measurements. 

We find that the scaling factors range from / sca i c ~ 
0.03 to 4.8, with a mean / SC aie,mean ~ 1-54 and a scat- 
ter on /scale, moan of Uf ~ 0.87. In «76% of the sources, 
the scaling factor is greater than unity, indicating that 
the variance is underestimated due to the sampling bias. 
Note that the Monte Carlo PSD is normalized by each 
source's variance, which is calculated using the source's 
light curve. Since the heavy binning could smear out 
high frequency variations, this method could result in an 
artificially smaller scaling factor. However, Allevato et 
al., in preparation, find a similarly small bias (/scale < 2) 
when sampling higher frequencies for a wide range of S/N 
ratios, gap lengths and sampling patterns. More impor- 
tantly, the scatter on the / sca ie factor calculated for each 
source is large (~ 40% — 190%), so individual measure- 
ments, even when corrected for bias, are likely to be poor 
estimates of the intrinsic variance. Therefore, variability 
properties of CDF-S galaxies are best considered in en- 
semble rather than on an individual basis. Bias-corrected 
excess variances (cr^ xs corr ) are listed in Table 2 and are 
used for all further analysis unless otherwise noted. 

4.2. Comparisons with XRB and ULX Variability 

The three most likely sources of X-ray variability are 
X-ray binary (XRB) populations, ultraluminous X-ray 
sources (ULXs), and accreting SMBHs. In this section, 
we show that the first two possibilities are not likely to 
dominate the measured galaxy variability. 

To examine the potential contribution of X-ray bi- 
nary populations to variability, we must first deter- 
mine the relative contributions of low-mass X-ray bi- 
naries (LMXB) and high-mass X-ray binaries (HMXB) 
to the hard (2-10 keV), galaxy- wide X-ray luminosity 
(-£'XRb)I3 Galaxy stellar mass (M*) scales the contribu- 
tion of (older) LMXBs, and star formation rate (SFR) 
scales the contribution of (you nger) HMXBs (e.g., Equa- 
tion 3 of ILehmer et all 1201 CD. M+ and S FR are calcu- 
lated for each galaxy in IXue et al.l ([2010) using the op- 
tical colors and the UV and IR lumin osities. Both the 
ILehmer et all (|2010D relations and the IXue et al.l (|2010f ) 
calculations ad opt the same in itial stellar mass fun ction 
(IKrou pa 2001). However, since ILehmer et all (|2010f) and 
IXue et all (|2010h use different formalisms for computing 
stellar masses (Bell et al. 2003 and Zibetti et al. 2009, 
respectively, which differ primarily due to their models 
of star formation history) , we a pply a corr e ction factor 
of 2.6 to the stellar masses from IXue etahl pOlCh . The 
SFR and M* values for each variable galaxy are listed in 
Table 2. 

By comparing the expected LMXB and HMXB contri- 
butions to the total luminosity, we can determine which 

14 We measure variability in the 0.5-8 keV band, but by limiting 
the comparison to the 2-10 keV luminosities, we limit the contri- 
bution of other sources of galaxy- wide X-ray emission (hot gas, 
supcrnovac, supernova remnants, and O-stars), which fade sharply 
above 2 keV and can be considered negligible. 



population ought to dominate the variability. In the vari- 
able galaxies, SFR ranges from 0.04 to 55 Mq yr _1 with 
a median of 2.2 Mq yr _1 ; M* ranges from 2 x 10 7 to 
6 x 10 n M Q with a median of 1.2 x 10 10 M o . We find that 
all but six of 20 variable galaxies are expected to have a 
larger HMXB contribution. Lhmxb(SFR)/L L mxb(-M*) 
ranges from 0.01 to 75.8 with a median of 2.6 (Ta- 
ble 2). HMXBs are ge nerally more variable than LMXBs 
(|Gilfanov et al.l 120041) . so unless LMXBs dominate the 
X-ray output of a galaxy, we neglect their contribution. 

To determine the variability expected from the 
HMXBs, we follow the relations in §4.2.3 of 
Gilfa nov et al.l (|2004f) . where the variability of the 
HMXB population is roughly determined by the 
galaxy's SFR. The following relations were obtained 
from Monte Carlo simulations with a power-law HMXB 
luminosity function with slope a = —1.6 and a cut-off 
luminosity at L cu t = 2xl0 40 



&rms,tot 
O~rms,0 



0.35 



+0.35 
-0.10 



ergs s : 
for SFR < 5 M Q yr" 



(4) 
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Q-30±g;J§ for SFR = 5-10 Mq yr^ 1 (5) 
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ffrm S ,0 SFR 1/2 



for SFR > 10 M Q yr" 
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Here, ay mSi o is the fractional rms (i.e., the square root of 
the excess variance) expected from an individual X-ray 
binary, which c an be as larg e as 20-30% on ~year 
timescales (e.g.. iGilfanovl [20101 ) . and a rms .tot is the to- 
tal variability. We take o~ rms q — 0.3 and calculate the 
upper limit on a rm s.tot- In the most extreme case, we 
find that the upper limit on variability expected from an 
HMXB population is <r| RB < 0.044. 

In six galaxies, LMXBs are expected to dominate the 
X-ray luminosity. Four of the six have log M* > 10.5 
Mq, and therefore have large enough stellar mass to fol- 

— K M S tellfr law (§ 4 - 4 - 2 of Gilfanov et al. 



low a <Trm " 

2004). For the rem aining two galaxies, w e follow the 
trend in Figure 12 of lGilfanov et all (|2004| ). The upper 
limit on variability expected from an LMXB population 
is CTxRB < 0.02. 

We find that XRB populations cannot explain the full 
extent of the X-ray variability for variable galaxies. Fig- 
ure [8] plots the distribution of o- 2 xs corr for CDF-S galax- 
ies, showing that all the variable galaxies exhibit vari- 
ability in excess of that expected from XRBs. The vari- 
able galaxies have a median c 2 XS;Corr /er^ RB w 42 (with- 
out the bias correction, o~^ s /a^ B ~ I 4 )' indicating that 
the contribution of an XRB population to the measured 
variability is small. 

We also consider whether one or more ultraluminous 
X-ray sources (ULXs) may dominate a galaxy's X-ray 
output. The nature of ULXs is still debated, but 
most likely involves accretion onto massive stellar black 
holes (30-100 -Mq); a few cases may involve accretion 
onto intermediate mass black holes (100-300 Mq) or 
beamed emission from 10-20 Mq black holes. Typical 
luminosities s pan Ln.s-skeV ~ 10 39 -10 41 ergs s _1 (e.g., 
ILehmer et aL1l2006l: ISwartz et al.l 120111 ). In a survey of 
1,441 X-ray point sources in 32 nearb y galaxies in the 
Chandra archive, iColbert et al.l (|2004l) found ULXs in 
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19 galaxies; the contribution of one or more ULXs to a 
galaxy's total X-ray point source luminosity ranged from 
7%-87%, with a median contribution of 43%. 

ULXs could potentially explain nine variable galaxies 
with £o.5_8kcV < 10 41 er S s s_1. However, since ULXs 
tend to be associa ted with star- forming regions (e.g., 
iSwartz et ail [20041 and occur more frequently in late- 
type /irregular gala xies than in early- type galaxies (e.g., 
I Walton et al.ll201lj) . the five variable galaxies with both 
£o.5- 8kcV < 10 41 ergs s _1 and late-type morphology (see 
£|5.2[) are more likely to host ULXs. 

One method of finding ULXs is to search for off- nuclear 
sources. We plot postage-stamp images (8" x 8") of 
the GOODS-S/GEMS HST U606-band for the variable 
galaxies in Figure HI The circle overplotted on each im- 
age has a radius 1.5 times the Chandra positional erro, 
w hich is calcula t ed at the 90% significance level. As 
in iLehmer et al.l ([2006?), X-ray sources offset from the 
galaxy nucleus by more than 1.5 times the positional er- 
ror are considered off-nuclear. 

We find one marginally off-nuclear source: 
033219.27-275406.7 (XID = 269). The primary 
optical source appears to be an early-type galaxy with 
blending toward the galaxy to the lower right, suggesting 
a merger. Both galaxies have similar redshifts (z = 0.960 
and 0.956, respectively); the first redshift is spectro- 
scopic and secur e (see $2]) and the second is photometric 
(|Xue et al.ll2010D . The high X-ray luminosity (£o.5-8kcV 
?» 2.2 xlO 42 ergs s -1 ) suggests that a ULX is likely not 
the dominant source of X-ray emission from this galaxy. 

One variable galaxy, 033230.00-274405.0 (XID = 418), 
was pre viously identified as being off-nuclear in the 1 Ms 
CDF-S (|Lehmer et al.l [20061. With the additional data 
from the 4 Ms CDF-S, the X-ray source position has 
been refined (with reduced uncertainty) to be consistent 
with the galaxy's nucleus. The off-nuclear source dis- 
cussed in the previous paragraph (XID = 269) was not 
detected in the 1 Ms CDF-S, so it was not considered by 
ILehmer et all (|2006f ). 

We cannot rule out that a ULX may lie too close to a 
galaxy's nucleus to be detected as an off-nuclear source. 
However, the possibility of ULXs in most variable galax- 
ies is mitigated by high X-ray luminosities and/or early- 
type morphology, so since we find no plausible off-nuclear 
ULXs, we conclude that ULXs are unlikely to dominate 
the emission from variable galaxies. Accretion onto a 
SMBH remains the best explanation of variable galaxies. 

5. SUPPORTING EVIDENCE FOR LLAGN 

We investigate the X-ray spectral shapes, the mor- 
phologies, and the optical spectral classifications of the 
variable galaxies for two purposes: (1) characterizing the 
variable galaxy population, and (2) determining whether 
their properties are consistent with those of LLAGN. 

5.1. X-ray Spectral Shape 

As discussed in fj2j Xue et all (|201 If ) calculate the ef- 
fective photon index for each source based on the ratio 
of count rates in the 2-8 keV and 0.5-2 keV bands (Ta- 
ble 2). For 10 variable galaxies detected in the soft band 
but not the hard band, lower limits are listed. For 4 low- 
count variable galaxies detected either in the full-band, 
soft-band, or both, no reliable effective photon index can 
be determined, so F e g is set to 1.40. 



Since 14 sources have poorly determined r e ff, we per- 
form a stacking analysis on all variable galaxies with 
< 150 net counts to determine an average photon index 
for the sample. The three highest count sources (with 
199.2, 275.6, and 608.9 net counts) are excluded from the 
stacking analysis since they could dominate the stacked 
signal, but the results remain the same within errors if 
these sources are included. 

Following the procedure in lLuo et al.l (|201 If ) . we calcu- 
late the soft and hard-band counts in a 3" diameter aper- 
ture for each source. The background is calculated by 
averaging the counts in 1,000 randomly placed, source- 
free apertures within a l'-radius circle around the source 
position. The individual source counts are summed and 
background is subtracted. Aperture corrections are ap- 
plied, averaged over all the observations weighted by ex- 
posure time (see Luo et al. 2011 for details), before cal- 
culating the band ratio. 

The stacked effective photon index for the 17 relevant 
variable galaxies is r sta ck ~ 1.93±0.13, which is consis- 
tent with the typical photon index for local Seyfert galax- 
ies (r ~ 1.8; e.g., Dadina et al. 2008) at the la level. 
Including all 20 variable galaxies, r sta ck ~ 1.82±0.08. 
While absorption may still be present in some individ- 
ual sources, the soft spectrum implied by the stacked ef- 
fective photon index indicates that the average variable 
galaxy is not heavily obscured. The X-ray luminosities 
are therefore intrinsically low and indicate that variable 
galaxies are most likely LLAGN. 

5.2. Galaxy Morphologies 

Postage-stamp images (8" x 8") of the GOODS- 
S/GEMS HST U606-band for CDF-S variable and non- 
variable galaxies are shown in Figures [9] and [TQl respec- 
tively, with Chandra error circles overlaid. Galaxies are 
classified by eye as late- type, early- type, irregular, or un- 
determined. Mergers are also visually classified based on 
blending and/or tidal tails between two or more galaxies. 

Since visual classifications are subjective and are par- 
ticularly unreliable for distant, poorly resolved galax- 
ies, we also appl y the color-magnitude relation given in 
I Bell et afl 1)20041 ) . where galaxies are considered part of 
the "red sequence" (i.e., early-type morphology) if they 
are redder than (M v - M v ) = -0.312 - 0.08M y - 0.51. 
Galaxies with bluer colors are considered part of the 
"blue cloud" (i.e., late-type morphology). Although the 
color-magnitude diagram leads to a more objective clas- 
sification scheme, there are nevertheless significant un- 
certainties, both in the rest-fra me magnitudes an d col- 
ors, and in the definition of the iBell et all (|2004| ) rela- 
tion. Moreover, the color classification is complicated 
by the fact that many of the galaxies lie in between the 
red sequence and blue cloud, in the so-called "green val- 
ley." In cases where the color and the visually classi- 
fied morphology disagree, we choose a final classification, 
preferring the visually classified morphologies in nearby, 
well-resolved galaxies, and the color classification in dis- 
tant and/or poorly resolved galaxies. Considering only 
galaxies classified as early or late-type, the visual and 
color classifications agree «54% of the time. Table 2 
lists each galaxy's morphological type as determined by 
visual classification, th e rest-frame magnitude and color 
ffr omlXue et al.ll2011l) . the color classification according 
to IBell et al.l (|2004l ). and the final classification. 



Based on the final classification from Table 2, we find 
that variability does not prefer one morphology type 
over the other. Variable galaxies have 40.0^g 1 g 3 % early- 
type morphology and 50. Ot^;^ late-type morphology, 
compared to 23. 6^3% and 51. 4^8%, respectively in 
non-variable galaxies (the errors are la binomi al er- 
rors). If we instead apply only the lBell et al.l (|2004l ) color 
classifications, then both the variable and non-variable 
galaxies prefer late-type morphologies (80.0lj^ 5 % and 

72.2±5 j%, respectively). The M* and SFR distributions 
for variable and non-variable galaxies show no significant 
difference (P KS > 10%). 

We find that the fraction of mergers among vari- 
able and non- variable galaxies is consistent within errors 
(9.1™% and 20.0±|;|%, respectively). 

5.3. Optical Spectroscopic Classification 

Most (18/20) var iable galaxies have o ptical spec- 
tral observations (ISzokolv et al.l l2004t iZheng et al.l 



20041: iMignoli et all 120051: iRavikumar et all l2007t 
Popesso et all 120091 : iSilverman et al.l l2010h ■ from which 



spectroscopic redshifts were determined. In all cases, 
the optical spectra are classified as galaxies, showing 
o nly narrow em is sion l ines or absorption lines. 

ISzokolv et all (|2004l) classify objects i n more detail 
Of the eight variable galaxies listed in the ISzokolv et all 
(2004) catalog, two have typical galaxy spectra showing 
only absorption lines. The remaining six are classified 
as having low-excitation emission lines consistent with 
H II region-type spectra. These objects would be classi- 
fied as normal galaxies based on the optical data alone 
as the presence of the AGN cannot be established. How- 
ever, one of these, 033222.78-275224.2 (XID = 312), has 
sufficient signal-to-noise in the optical spectrum to mea- 
sure line ratios. Thi s object is classified as a LINER by 
ISzokolv et all (l2004h via the line ratio diagnostics given 
bv lHo fc Sargent! (|1993D . 

6. GALAXY VS. AGN VARIABILITY 

6.1. The Variability- Luminosity Anti- Correlation 

In this section, we investigate how galaxy variability 
compares to AGN variability. We first confirm a signifi- 
cant anti-correlation between excess variance and X-ray 
luminosity among AGN, as seen in previous work (e.g. , 
Barr fc Mushotzkyj 119861: lLawrence fc Pa padaki s 119931 : 
Nandra et all Il997t iHawkinsI 120001 : iPaolillo et alH2004h . 
The anti-correlation is plotted in Figure [TT] for signifi- 
cantly variable sources as a^ xs corr vs. Lo.5-8kcV- The 
exce ss variance includes the bias correction discussed in 
N4. 11 The rest-frame X-ray luminosities are calculated 
differently for AGN and galaxies, as described in Sj2j The 
Spearman rank correlation coefficient (p s = —0.44) shows 
the correlation is significant at P s w 10 -9 (5.9<r) for AGN 
only (red circles). The correlation increases in signifi- 
cance to P s « 8 x 10~ 12 (6.4cr) if all variable sources 
above Lo.5-8keV = 10 41 ergs s _1 are considered, includ- 
ing those classified as galaxies (black circles and stars). 
No significant correlation is found if the sample is limited 
to only the variable galaxies, which are discussed further 
in £16.21 Note that the correlation coefficient and best-fit 
line parameters (given below) remain consistent within 



errors if calculated using the excess variance uncorrected 
for sampling bias (ct 2 xs ). 

Since both variability and luminosity may depend on 
other parameters, such as black hole mass and/or accre- 
tion rate, we fit the data using SIXLIN .PRO, an IDL pro- 
gram adapted from llsobe et al.l (|1990l ) . A least-squares 
bisector fit to sources with Lo.5-8 kcV > 10 41 ergs s _1 , 
weighted by the uncertainties in cr 2 xs corr , results in 

a 2 nxStCotI = (25.8±2.6) L ? ° 5 ^ v 06 . The slope and in- 
tercept are consistent within errors if the weights are 
not included. Large blue squares in Figure [TT] show the 
weighted means for each luminosity bin; as expected, 
these are consistent with the weighted least-squares bi- 
sector. Note that since weighted means will most closely 
follow data with small error bars, they are weighted heav- 
ily downward in this log-log plot. Limiting the sample 
further to sources with Lq 5-8keV > 10 42 ergs s _1 results 
in ^ xs ,corr = (31.9±2.7) L"° 5j !™_Theslopes are both 
consistent with the results of lNandra et all (|1997l ) , where 
a weighted least-squares bisector fit results in er 2 xs cx 
^a-iotev ' 03 (Nandra et al. 1997 do not give the normal- 
izat ion of their relation) . This result is notable because 
the iNandra et all (|1997f) data sampled shorter (hour- 
day) timescales compared to the months-years timescales 
sampled by the CDF-S. 

The slope presented in this paper is significantly flat- 
ter than that found in Paolillo et al. (2004; cj 2 xs 
oc gi^y 23 ), which included non- varying sources. 

We choose not to include upper limits for non-varying 
sources in our analysis because the assumptions under- 
lying survival analysis, which have been successfully ap- 
plied to deal with censored data in other astronomical 
situations, do not apply here because: (1) the excess 
variance measurements of most sources lie near the de- 
tection limit, (2) a large percentage of sources do not 
have detected variability («50% of AGN and «78% of 
galaxies), and (3) a significant percentage of sources 
with no detected variability, especially those classified 
as galaxies, may truly not be variable (i.e., <r 2 xs ~ 0). 
By not including censored data in this paper's measure- 
ments, we likely bias the measured slopes and possibly 
the significance of the anti-correlation. Nevertheless, the 
variability-luminosity anti-correlation has been observed 
to follow a model based on SMB H mass and accretion 
rate (e.g.. iPapadakis et al.l l2008bh . suggesting that the 
anti-correlation is real, though probably not linear. The 
model is discussed further in H6.3I 

We check for other potential biases that may affect the 
luminosity-variability anti-correlation. The log-log plot 
of Figure QTJ has the disadvantage of "hiding" sources 
with negative tr 2 xs corr values. Sources with low flux val- 
ues, and hence higher scatter in tr 2 xs corr will therefore 
appear to have stronger variability, since the <7 2 xs corr 



values scattered to negative values will be hidden. To 
check for this bias, we remove high scatter sources with 
err(cr 2 xs corl .) > 0.1, excluding all but 38 sources, and 

find that the er 2 xs corr — £o.5-8keV correlation remains sig- 
nificant at P s 0.9% (2.6ct). The best-fit slope and 
intercept remain consistent at the la level. 

The flux-limited nature of the CDF-S survey presents 
another potential bias. Since luminosity correlates with 
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redshift, and intrinsic variability timescales decrease with 
redshift, the decrease in variability strength as luminos- 
ity increases could, in principle, simply reflect the fact 
that shorter timescales are studied at higher redshifts, 
and therefore exhibit less variability due to the red-noise 
nature of AGN light curves. To check for a possible red- 
shift bias, we examine the a^ xs corr -Lo.5-8keV correlation 
using a sub-sample within a narrow redshift range (0.55 
< z < 0.75). This redshift range selects 10% of the to- 
tal sample, covers luminosities from Lo.5-8kcV ~ 10 415 
to 10 44 ergs s _1 , and results in minimal differences in 
rest-frame timescales. The Spearman rank correlation 
coefficient for this sub-sample remains significant at P s 
« 0.8% (2.6cr); the slope and intercept are consistent at 
the lcr level. 

The above test also addresses another potential bias 
due to the redshift range of the sample. The a„ xs val- 
ues listed in Table 2 measure the variability strength in 
the observed frame, so they will sample different energy 
bands depending on the source's redshift (see ^4.11 for 
discussion). This could introduce bias if the variablity 
amplitude changes with increasing energy: from 0.5 — 8 
keV atz = 0tol — 16 keV at z = 1. The slope and inter- 
cept remain the same within the narrow redshift range 
tested above, suggesting that any such bias does not have 
a significant effect on the sample. 

6.2. Suppressed Variability in LLAGN 

Variable sources with luminosities less than Lo.5_8keV 
= 10 41 ergs s -1 tend to fall significantly below the ex- 
trapolated linear relation by factors of «6-80 (median 
factor of ~24), indicating a drop in variability relative to 
the linear relation on long timescales for LLAGN. This 
"suppressed" variability can be shown to be intrinsic to 
properties of AGN variability rather than due to dilution 
by unrelated XRB populations. 

In we showed that galaxy variability cannot be 

attributed solely to XRB populations. We now check 
whether the XRB contribution could nevertheless dilute 
the observed variability by estimating how much XRBs 
are expected to contribute to the total hard (2-10 keV), 
galaxy- wide X-ray luminosity (Lxrb; see ^4.21 Lehmer 
et al. 2010). 

We compare Lxrb to the measured, intrinsic X-ray 
luminosity. First, Lo.5-8keV is converted to L2-iokeV! 
using the i ntrins ic photon index of T = 1.8 adopted by 
IXue et al.l ([201 ID for AGN, and the stacked photon in- 
dex r sta ck ~ 1.9 for galaxies (see Figure IT21 plots 
the measured 2-10 keV luminosity against that expected 
from an XRB population for variable galaxies (red stars) 
and non- variable galaxies (green squares). For reference, 
we plot in the same figure 36 9 CDF-S AGN (ope n or- 
ange circles), 32 local galaxies ([Colbert et al.ll2004l ). and 
20 local luminous infrare d galaxies, which ar e likely to 
be actively star- forming ([Lehmer et al.l [2010). The red 
solid line shows unity, while the dotted li nes mark the 
dispersion observed in lLehmer et al.l (|2010f ). 

A K-S test shows no significant difference between the 
relative XRB contribution in variable and non-variable 
galaxies. The CDF-S galaxies on average tend to lie 
above unity, with a median i2_iokev/-^XRB ~ 4.6. The 
X-ray excess in the non- variable galaxy population is per- 
haps not surprising given the flux-limited nature of the 



CDF-S survey, which is more likely to detect the objects 
at the bright end of the galaxy X-ray luminosity function, 
and may result in a high percentage of "contamination" 
by AGN. 

The median X-ray excess for variable galaxies 
(L 2 -iokcV,tot/£2-iokcV,XRB ~ 9 - 2 ) suggests that XRBs 
contribute ~11% of the 2-10 keV luminosity for the aver- 
age variable galaxy. For six variable galaxies, the total X- 
ray emission is consistent with tha t expected from XRB 
emission within the scatter of the ILehmer et al.l ([20100 
relation (see Fig. [T2]), suggesting that dilution may be 
more significant in these sources. Three of these sources 
have variability consistent with the linear variability- 
luminosity relation, while three have suppressed variabil- 
ity (filled black stars in Fig. [TT1) . Dilution by XRB vari- 
ability may therefore play a role in suppressed variability, 
but cannot fully explain the extent to which er^ xs corr is 
suppressed at low luminosities. 

An alternative possible explanation for the suppressed 
variability at low luminosit ies is a change in accretion 
structure. iPtak et all (|1998ft found a similar drop in vari- 
ability strength below L2-iokcV ~ 2xl0 41 ergs s _1 in a 
sample of LLAGN and LINERs observed with ASCA on 
variability timescales of less than a day. The authors 
hypothesized that a radiatively inefficient accretion flow 
(RIAF, e.g., Yuan & Narayan 2004) could be responsi- 
ble for suppressed short-timescale variability at low lu- 
minosities due to the larger extent of the X-ray source. 
This scenario would not obviously explain the reduced 
variability on ~month-year timescales seen here. RIAF 
models also predict a hard X-ray photon index due to 
the lack of an optically thick accretion disk, which pro- 
vides the soft X-ray photons. The stacked X-ray photon 
index for variable galaxies (r sta ck ~ 1.93±0.13; i )5.ip is 
inconsistent with this predict io n. 

Studies since IPtak et all (|1998f ) have found evi- 
dence both for "s uppre ss ed" variability in LLAGN 
(e.g.. IPtak et al.l [2001 IMarkowitz fe Uttlevl [2001 
Papadakis et a l. 2008a) and against it (e.g., IBinder et all 
20091: iPian et all [2010). Similarly, objects such as 
narrow-line Seyfcrt 1 (NLS1) galaxies and the dwarf 
Seyfert NGC 4395 {M BH « 3.6 x lO 5 ^) exhibit 
"exce s s" vairability for th eir luminosity ([Boiler et all 
119961 : llwasawa et al.l I2O1O0 . However, when plotting 
variability against mass instead of lumin osity (e.g., 
IPapadakis et al.l I2008at IMiniutti et "aLll2009D . such dis- 
crepancies disappear, with residual diff erences possibly 
due to varying accretion rate s (e.g., [McHard v et all 
[2001 [Markowitz fc Uttlevl [20051) . 

To investigate the role of SMBH mass {Mbh) and ac- 
cretion rate (normalized by the Eddington rate; rh = 
M / 'M Edd) , we have obtained rough estimates for all vari- 
able sources. Masses and Eddington ratios for all galax- 
ies, variable and non-variable, are listed in Table 2. The 
SMBH masses are estimated via the scaling rel ation be- 
tween Mbh and absolute if-band magnitude (|Grahaml 
[2001 : 

log— = -0.37(±0.04)(M K + 24) + 8.29(±0.08) (7) 

The total absolute rest-frame if - band mag n itudes 
are derived from SED fitting in IXue et all (|2010fi 
with a random scatter of <0.3 mag. An X-ray 
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luminosity-dependent correction factor (Equation 1 of 
iVasudevan et al.|[2009l ) corrects for nuclear emission. We 
assume that the host galaxy is b ulge-dominated, a valid 
assumption for mo st AGN (e.e.. iKaufhnan et alJl2003t 
IGrogin et al.1 [2005). Several variable galaxies, how- 
ever, are not bulge-dominated (nine are late-type; see 
Table 2), so their black hole masses may be overes- 
timated. We apply a luminosity-dependent bolomet- 
ric correct ion (/t2-inkRv) to est imate the bolometric lu- 
minosity (Marconi et al.l l2004[ ) and calculate the Ed- 
dington ratio (L bol /L Edd = [K 2 -iokoV ^2 -iok e y]/[l-25 x 
10 38 (Mrh/Mc?))]). Note, that the iMarconi et"all 
(|2004l ) correction is calculated explicitly for L 2 -iokeV > 
10 42 ergs s _1 ; we extrapolate this relation down to the 
lower luminosities of the variable galaxy sample. 

Both the Mbh and Lboi estimation techniques are 
known to have large dispersions. The Mbh-Mk scaling 
relation has a total scatter of 0.33 dex, and additional un- 
certainty will come from the luminosity-dependent cor- 
rection fo^jruclea£_ejnissi£n, which is based on template 
SEDs ijVasudevan et al.1l2009D . In addition, the assump- 
tion that all the variable galaxies are bulge-dominated 
will produce additional scatterJ3 The bolometric cor- 
rection, too, has large scatter due to the intrinsic dis- 
persion in the SED. The uncertainty in Lboi due to SEP 
dispersion is ~20% f or luminous AGN (jElvis et al.lll994t 
iRichards et aL|[2006h . There is some debate regarding 
the similarity between the SEDs of LLAGN and lumi- 
nous AGN (e.g., gU Ho 1999 and Ho 2002 vs. Maoz et 
al. 2007); however, the luminous AGN SED dispersion 
likely serves as a lower limit to dispersion in LLAGN 
SEDs. 

Variable galaxies tend to have lower accretion rates 
({rh) = 4 x 10" 4 ) and masses ({Mbh) = 2.6 x 1O 7 M ) 
than variable AGN ({rh) = 9 x 10~ 3 ; (M BH ) = 6.2 x 
10 7 M Q ), where we have limited the AGN sample to z < 1 
for purposes of comparison. A K-S test shows that the 
differences in the m and Mbh distributions are signifi- 
cant: P K s ~ 1.7xl0~ 5 and P^s « 0.002, respectively. 
The properties of variable AGN and galaxie s are consis- 
tent w ith the range of estimates made by iBabic et al.1 
(2007) for X-ray selected, z < 1 AGN in the 1 Ms 
CDF-S, which span rh ~ 10" 5 - 1 (median » 0.001) and 
M BH ~ 10 6 - 10 10 M Q (median » 10 8 M w ). 

Unlike prev ious stu dies (lO'Neill et af] 120051: 

iPapadakis et al.l I2008at iMiniutti et al.1 12009ft . we 
find no significant (anti-)correlation between a^ xs and 
Mbh, most likely due to the large scatter in a^ xs and 
Mbh measurements, combined with the narrow range 
of masses probed (since most sources lie between Mbh 
~ 10 7 -10 9 Mq). It is not surprising, therefore, that we 
find a significant anti-correlation between a^ xs and rh 
(Pks ~ 1.6 x 10 -11 , 6.3ct), plotted in Figure fT3] which is 
likely an artefact of the (7„ xs -Lo.5-8keV anti-correlation. 
Nevertheless, it is interesting that in Figure [T3l the vari- 
able galaxies (black circles) connect smoothly with the 
more luminous AGN (red circles), with no discrepancy 
in variability strength. This suggests that the factor 
of 22.5 difference in rh may explain the "suppressed" 

15 The scaling relation in Kormcndy & Gcbhardt (2001), for ex- 
ample, has a much larger scatter of 0.56 dex largely because of 
poor bulge/disc separation (Graham 2007). 



variability in variable galaxies. 

6.3. Comparing the Variability-Luminosity Relation 
with Empirical Models 

A number of re c ent studies (e. g.JMcHardv et al1l2004l: 
Papadakisl [200l lO'Neill et all 120051: IPapadakis et all 
2008b) have shown that X-ray variability may be de- 
termined by a combination of Mbh and rh, explain- 
ing the observed variability-luminosity relation. AGN 
light curves appear to be universally described by a 
broken power-law PSD function, where the break fre- 
quency depends on mass and accretion rate: Vbf — 
0.029r]m(MBH /10 6 M Q ) -1 , wher e r? is the accretion ef - 
ficiency, assumed to be rj — 0.1 (McHa rdv et alj 12006). 
The excess variance is equivalent to the integral of the 
PSD between the minimum and maximum frequencies 
sampled by the data, so as long as the break frequency 
falls within this range, the excess variance at a given lu- 
minosity will decrease with increasing Mbh and increase 
with increasing rh. 

We compare the Mbh and rh values estimated in %.2\ 
to the variability measured for the variable AGN and 
galaxy populations. Assuming a universal PSD func- 
tion, we derive the variability-luminosity relations ex- 
pected for the CDF-S sampling pattern, given a range of 
SMBH masses and accretion rates fe.g..|Papadaki s~et al.1 
2008b, and references therein). The bolometric luminos- 
ity, which is calculated directly from Mbh and rh, is con- 
yerted to X-ray lum inosity via a bolometric correction 
(M arconi et al.l 12004). Where the break frequency lies 
outside the timescales sampled by the data (depending 
on the combination of rh and Mbh), the excess variance 
will remain constant. 

The model variability-luminosity relations are plotted 
in Figure [14] for the average accretion rates for variable 
(rh = 9 x 10~ 3 ; dashed line) and non-variable (rh — 
4 x 10 -4 ; dash-dotted line) galaxies covering Mbh = 
10 5 -10 1 M (low to high L .5_ 8 kcv)- Both variable 
populations are limited to z < 1 for comparison, and 
the relations are calculated at z = 0.5 and z = 0.7, the 
median redshifts for variable galaxies and AGN, respec- 
tively; using model redshifts z = or 1 resulted in neg- 
ligible changes. The shape of the observed variability- 
luminosity anti-correlation, including the plateau at low 
luminosity, is roughly reproduced by the model. The 
plateau occurs when the break frequency lies outside 
the range of timescales sampled by the data. Unfortu- 
nately, the long timescales (especially the long minimum 
timescale of 4 Ms; see t j4.1[) and the large scatter prevent 
the models from distinguishing between the significantly 
different accretion rates estimated for variable galaxies 
and AGN. 

Note that most sources are more variable than pre- 
dicted by the model. We note some possible sources of 
bias: (1) The normalization of the PSD function used 
to calculate t he models is ba sed on a small sample of 
nearby AGN (jPapadakis 2003) and may therefore not be 
representative of CDF-S AGN out to z « 1. (2) Both 
the models and the bias correction applied to the mea- 
sured excess variance (t j4.ip depend on a universal bro- 
ken power-law PSD function, but the slopes of the PSD 
may vary between individual sources, and some sources 
may even have a second break at shorter frequencies 
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fe.g.. lMcHardv et al.ll2007h . (3) Peculiar sources, such as 
NLSls, exhibit higher variability^, so the variable sam- 
ple may suffer from selection effects due to the higher 
sensitivity to more variable sources ( H4.ll) . 

7. CONCLUSIONS AND FUTURE WORK 

We have investigated the variability properties of 
sources classified as galaxies in the 4 Ms CDF-S by di- 
viding the observations into ~1 Ms epochs. We find the 
following main results: 

1 . Using a x 2 -based Monte Carlo simulation appropri- 
ate for low-count sources, 20 of 92 galaxies are selected 
as variable, as well as 185 of 369 AGN. Variability is ef- 
fective in selecting AGN that do not meet other AGN 
selection criteria, such as cuts based on X-ray luminos- 
ity, X-ray-to-optical flux ratio, or excess X-ray emission 
compared to radio emission expected from star-forming 
galaxies (see $3]). If all 20 variable galaxies are accreting 
SMBH, t hen the common ly used selection criteria em- 
ployed bv lXue et al.l (l20"H miss 20 of 389 AGN («5.1%), 
and 19 of 54 LLAGN («35.2%; where LLAGN are arbi- 
trarily defined as having £o.5-8kcV < 10 42 ergs s _1 ). Due 
to the low counts measured for most galaxies, it is likely 
that some non-variable galaxies may still host an AGN 
but not exhibit detectable variability, so these fractions 
serve as lower limits. 

2. We measure variability strength with the normal- 
ized excess variance, accounting for measurement error, 
red-noise scatter, and bias due to the CDF-S sampling 
pattern. Comparing the bias-corrected excess variance 
to that expected from XRB populations, we find that 
XRBs cannot explain galaxy variability. 

3. The possibility of energetically significant ULXs in 
variable galaxies is mitigated by high X-ray luminosities 
and/or early- type morphology. In addition, we find no 
plausible off- nuclear ULXs. 

4. Galaxy X-ray variability is most likely associ- 
ated with accretion onto a relatively unobscured SMBH. 
Though some absorption may be present, the individ- 
ual and stacked photon indices show no indication of 
heavy obscuration and are consistent with the typical 
AGN spectral shape (r sta ck ~ 1.93±0.13). 

5. We confirm a significant anti-correlation between 
excess variance and X-ray luminosity, and find the slope 
and intercept to be consistent with what is reported in 
the literature for shorter timescales. We show that the 
sampling bias induced by the CDF-S observing pattern 
does not affect the anti-correlation slope or significance. 

6. Low-luminosity sources (io.5-8kcV < 10 ergs s _1 , 
largely variable galaxies) have "suppressed" variability 
compared to the extrapolated linear relation between ex- 
cess variance and X-ray luminosity. This may be ex- 
plained by their lower accretion rates. Variable galax- 
ies sample a different mass-accretion rate space than the 
rest of the AGN population, with significantly lower ac- 
cretion rates ((m) = 4 x 10~ 4 ) and masses {(Mbh) = 
2.6 x 1O 7 M ) than variable AGN ((m) = 9 x 10~ 3 ; 
(Mbh) = 6.2 x 10 7 M Q ) in the same redshift range 
(z<l). 

16 A notable exception is the NLS1 galaxy Akn 564, which is 
the only AGN observed to date to have a second, low-frequency 
break in it s power spectrum, lea ding to lower variability on long 
timescales HMcHardv et al.ll2u07n . 



7. We find that an empirical model based on a 
universal broken power-law PSD function, where the 
break frequency depends on SMBH mass and accretion 
rate, roughly reproduces the shape of the variability- 
luminosity anti-correlation. However, the normalization 
of the model is low compared to the data, suggesting sys- 
tematic bias either within the data (e.g., selection effects) 
or the model (e.g., incorrect normalization for LLAGN 
at z < 1). 

We have shown that the variability measured by deep 
X-ray surveys is an effective technique for selecting cos- 
mologically distant LLAGN. Extending the CDF-S to 
even longer exposures would enable detection of variabil- 
ity in both fainter and less variable sources, thus allowing 
better characterization of the properties and abundance 
of the LLAGN population. Within the 4 Ms CDF-S, 
smaller time bins could be used to search for variability 
and better characterize the variability presented here. 

In addition, follow-up optical spectroscopy is necessary 
to ch aracterize the LLAGN population presented here. 
While ISzokolvet al.l ((2004) were able to classify one ob- 
ject as a LINER, the spectra for the other sources did 
not have sufficient signal-to-noise to measure line ratios, 
so these objects would require either deeper optical spec- 
troscopy or stacking of the currently available spectra. 

We thank P. Uttley, E. Feigelson, and C. Saez for help- 
ful discussions on variability statistics. We also thank M. 
Gilfanov for helpful discussions regarding XRB variabil- 
ity. We thank the referee for constructive comments. 
We acknowledge the financial support of NASA ADP 
grant NNX10AC99G (M.Y., W.N.B., Y.Q.X.) and Chan- 
dra X-ray Center grant SP1-12007A (W.N.B., Y.Q.X.). 
M.P. acknowledges support from PRIN-2009 by the Ital- 
ian MIUR. F.E.B. acknowledges support from Programa 
de Financiamiento Basal and CONICYT-Chile under 
gra nts FO NDECYT 1101024, alma-conicyrSTT OOOOl, 
and fondap-cata 15010003. B.D.L. acknowledges finan- 
cial support provided by the Einstein Fellowship Pro- 
gram. D.M.A. thanks the Science and Technology Facil- 
ities Council (STFC) for support. 

REFERENCES 

Alexander, D. M., Bauer, F. E., Chapman, S. O, Smail, I., Blain, 

A.W., Brandt, W. N., & Ivison, R. J. 2005, ApJ, 632, 736 
Alexander, D. M. et al. 2008, ApJ, 687, 835 
Allcvato, V. et al. 2011, in preparation 
Almaini, O. et al. 2000, MNRAS, 315, 325 

Babic, A., Miller, L., Jarvis, M. J., Turner, T. J., Alexander, D. 

M. & Croom, S. M. 2007, A&A, 474, 755 
Balestra, I., et al. 2010, A&A, 512, 12 
Barr, P. & Mushotzky, R. F. 1986, Nature, 320, 421 
Barmby P. et al. 2006, ApJ, 650, L45 

Barth, A. J., Ho, L. O, Filippenko, A. V., & Sargent, W. L. W. 

1998, ApJ, 496, 133 
Bauer, F. E., Alexander, D. M., Brandt, W. N, Schneider, D. P., 

Treister, E., Hornschemeier, A. E., & Garmire, G. P. 2004, AJ, 

128, 2048 

Bell, E. F., Mcintosh, D. H., Katz, N., & Weinberg, M. D. 2003, 
ApJS, 149,289 

Bell, E. F., Wolf, C, Meisenheimer, K., Rix, H.-W., Borch, A., 
Dye, S., Kleinheinrich, M., Wisotzki, L.& Mcintosh, D. H. 
2004, ApJ, 608, 752 

Binder, B., Markowitz, A. & Rothschild, R. E. 2009, ApJ, 691, 
431 

Blundell, K. M., Rawlings, S., & Willott, C. J. 1999, AJ, 117, 677 



12 



Boiler, T., Brandt, W. N. & Fink, H. 1996, A&A, 305, 53 
Boutsia, K., Leibundgut, B., Trevese, D. & Vagnetti, F. 2009, 
A&A, 497,81 

Brandt, W. N. & Hasinger, G. 2005, ARA&A, 43, 827 

Broos, P. S., Townsley, L. K., Feigelson, E. D., Getman, K. V., 

Bauer, F. E., & Garmire, G. P. 2010, ApJ, 714, 1582 
Caldwell, J. A. R., et al. 2008, ApJS, 174, 136 
Cameron, E. 2011, PASA, 28, 128 
Cardamone, C. N., et al. 2010, ApJS, 189, 270 
Colbert, E. J. M., Heckman, T. M., Ptak, A. F., Strickland, D. K. 

& Weaver, K. A. 2004, ApJ, 602, 231 
Cooper, M. C. et al. 2012, MNRAS, submitted, arXiv:1112.0312 
Dadina, M. 2008, A&A, 485, 417 
Damen, M., et al. 2011, ApJ, 727, 1 
Elvis, M. et al. 1994, ApJS, 95, 1 

Ferrarese, L. & Ford, H. 2005, Space Sci. Rev., 116, 523 
Freeman, P. E., Kashyap, V., Rosner, R., & Lamb, D. Q. 2002, 

ApJS, 138, 185 
Gallo, L. C, Boiler, Th., Brandt, W. N., Fabian, A. C. & 

Vaughan, S. 2004, MNRAS, 355, 330 
Gallo, L. C, Brandt, W. N., Costantini, E. & Fabian, A. C. 2007, 

MNRAS, 377, 1375 
Gehrels, N. 1986, ApJ, 303, 336 
Giavalisco, M., et al. 2004, ApJ, 600, L93 

Gilfanov, M., Grimm, H. J. & Sunyaev, R. 2004, MNRAS, 351, 
1365 

Gilfanov, M. 2010, LNP, 794, 17 

Goncalves, T. S., Steidel, C. C, & Pettini, M. 2008, ApJ, 676, 816 
Gonzalez-Martin, O., Papadakis, I., Reig, P. & Zezas, A. 2011, 

A&A, 526, 132 
Graham A. 2007, MNRAS, 379, 711 
Grazian, A., et al. 2006, A&A, 449, 951 
Grogin N. A. et al. 2005, ApJ, 627, 97 
Gultekin, K., et al. 2009, ApJ, 698, 198 
Hawkins, M. R. S. 2000, A&A, 143, 465 
Hickox, R. C. & Markevitch, M. 2006, ApJ, 645, 95 
Hopkins, P. F. & Hernquist, L. 2009, ApJ, 698, 1550 
Hornschemeier, A. E. et al. 2003, ApJ, 126, 575 
Ho, L. C, Filippenko, A. V., & Sargent, W. L. W. 1993, ApJ, 

417, 63 

Ho, L. C, Filippenko, A. V., & Sargent, W. L. W. 1997, ApJS, 

112, 315 
Ho L. C. 1999, ApJ, 516, 672 
Ho L. C. 2002, ApJ, 564, 120 
Houck, J. R. et al. 2005, ApJL, 622, 105 

Isobe, T., Feigelson, E. D., Akritas, M. G. & Babu, G. J. 1990, 

ApJ, 364, 104 
Iwasawa, K., Tanaka, Y. & Gallo, L. C. A&A, 514, 58 
Jakobsen, P., Jansen, R. A., Wagner, S., & Reimers, D. 2003, 

A&A, 397,891 
Kauffmann, C, et al. 2003, MNRAS, 346, 1055 
Kellermann, K. I., Sramek, R., Schmidt, M., Shaffer, D. B. & 

Green, R. 1989, AJ, 98, 1195 
Kormendy J. & Gebhardt K., 2001, in Martel H., Wheeler J. C, 

eds, The 20th Texas Symposium on Relativistic Astrophysics. 

Am. Inst. Phys., New York, 586, 363 
Kroupa, P. 2001, MNRAS, 322, 231 
Komatsu, E., et al. 2011, ApJS, 192, 18 
Lawrence, A. & Papadakis, I. 1993, ApJ, 414, 85L 
Le F'evre, O., et al. 2004, A&A, 428, 1043 

Lehmer, B. D., Brandt, W. N., Hornschemeier, A. E., Alexander, 
D. M., Bauer, F. E., Koekemoer, A. M., Schneider, D. P. & 
Steffen, A. T. 2006, AJ, 131, 2394 

Lehmer, B. D., et al. 2008, ApJ, 681, 1163 

Lehmer, B. D., Alexander, D. M., Bauer, F. E., Brandt, W. N., 
Goulding, A. D., Jenkins, L. P., Ptak, A. & Roberts, T. P. 
2010, ApJ, 724, 559 

Luo, B., et al. 2008, ApJS, 179, 19 

Luo, B., et al. 2010, ApJS, 187, 560 

Luo, B., et al. 2011, ApJ, 740, 37 

Marconi, M., Risaliti, G., Gilli, R., Hunt, L. K., Maiolino, R. & 

Salvati, M. 2004, MNRAS, 351, 169 
Markowitz A. & Uttley P. 2005, ApJ, 625, L39 
McHardy, I. M., Papadakis, I. E., Uttley, P., Page, M. J. & 

Mason, K. O. 2004, MNRAS, 348, 783 
McHardy, I. M., Koerding, E., Knigge, C, Uttley, P. & Fender, 

R. P. 2006, Nature, 444, 730 



McHardy, I. M., Arevalo, P., Uttley, P., Papadakis, I. E., 

Summons, D. P., Brinkmann, W. & Page, M. J. 2007, MNRAS, 
382, 985 

Maoz, D. 2007, MNRAS, 377, 1696 

Maoz, D., Nagar, N. M., Falcke, H. & Wilson, A. S. 2005, ApJ, 
625, 699 

Maoz, D., Filippenko,A.V., Ho, L.C., Rix, H.-W., Bahcall, J. N., 
Schneider,D. P., & Macchetto, F. D. 1995, ApJ, 440, 91 

Martinez-Sansigre, A., Rawlings, S., Lacy, M., Fadda, D., 

Marleau, F., Simpson, C, Willott, C. J., & Jarvis, M. J. 2005, 
Nature, 436, 666 

Mignoli, M., et al. 2005, A&A, 437, 883 

Miller, N. A., Fomalont, E. B., Kellermann, K. I., Mainieri, V., 
Norman, C, Padovani, P., Rosati, P., & Tozzi, P. 2008, ApJS, 
179, 114 

Miniutti, G., Ponti, G., Greene, J. E., Ho, L. C, Fabian, A. C. & 

Iwasawa, K. 2009, MNRAS, 394, 443 
Morokuma, T. 2008, ApJ, 676, 121 

Nandra, K., George, I. M., Mushotzky, R. F., Turner, T. J. & 

Yaqoob, T. 1997, ApJ, 476, 70 
O'Neill, Paul M., Nandra, K., Papadakis, I. E. & Turner, T. J. 

2005, MNRAS, 358, 1405 
Paolillo, M., Schreier, E. J., Giacconi, R., Koekemoer, A. M., & 

Grogin, N. A. 2004, ApJ, 611, 93 
Papadakis I. E. 2004, MNRAS, 348, 207 

Papadakis, I. E., Ioannou, Z., Brinkmann, W. & Xilouris, E. M. 

2008a, A&A, 490, 995 
Papadakis, I. E., Chatzopoulos, E., Athanasiadis, D., Markowitz, 

A. & Georgantopoulos, I. 2008b, A&A, 487, 475 
Peterson, K. C, Gallagher, S. C, Hornschemeier, A. E., Muno, 

M. P. & Bullard, E. C. 2006, AJ, 131, 133 
Pian, E., Romano, P., Maoz, D., Cucchiara, A., Pagani, C. & La 

Parola, V. 2010, MNRAS, 401, 677 
Popesso, P., et al. 2009, A&A, 494, 443 

Ptak, A., Yaqoob, T., Mushotzky, R., Serlemitsos, P. & Griffiths, 

R. 1998, ApJ, 501, 37L 
Ptak, A., Terashima, Y., Ho, L. C. & Quataert, E. 2004, ApJ, 

606, 173 

Rafferty, D. A., Brandt, W. N., Alexander, D. M., Xue, Y. Q., 
Bauer, F. E., Lehmer, B. D., Luo, B., & Papovich, C. 2011, 
ApJ, 742, 3 

Ravikumar, C. D., et al. 2007, A&A, 465, 1099 
Richards, G.T. et al. 2006, ApJS, 166, 470 

Sarajcdini, V., Koo, D. C, Klesman, A. J., Laird, E. S., Perez 

Gonzalez, P. G. & Mozena, M. 2011, ApJ, 731, 97 
Scheuer, P. A. G. 1995, MNRAS, 277, 331 

Shankar, F., Salucci, P., Granato, G. L., De Zotti, G., & Danese, 

L. 2004, MNRAS, 354, 1020 
Silverman, J. D., et al. 2010, ApJS, 191, 124 

Supper R., Hasinger G., Pietsch W., Trmper J. Jain A., Magnier 
E.A., Lewin W.H.G. & Van Paradijs J. 1997, A&A, 317, 328 

Swartz, D. A., Ghosh, K. K., Tennant, A. F., & Wu, K. 2004, 
ApJS, 154, 519 

Swartz, D. A., Soria, R., Tennant, A. F. & Yukita, M. 2011, ApJ, 
741, 49 

Szokoly, G. P., et al. 2004, ApJS, 155, 271 
Taylor, E. N., et al. 2009, ApJS, 183, 295 
Timmer J. & Koenig M. 1995, A&A, 300, 707 
Tozzi, P., et al. 2001, ApJ, 562, 42 
Treister, E., et al. 2009, ApJ, 693, 1713 

Trevese, D., Boutsia, K., Vagnetti, F., Cappellaro, E. & Puccetti, 

S. 2008, A&A, 488, 73 
Vaughan, S., Boiler, T., Fabian, A. C, Ballantyne, D. R., Brandt, 

W. N., & Trumper, J. 2002, MNRAS, 337, 247 
Vaughan, S., Edelson, R., Warwick, R. S. & Uttley, P. 2003, 

MNRAS, 345, 1271 
Vaughan, S., Fabian, A. C, Ballantyne, D. R., De Rosa, A., Piro, 

L. & Matt, G. 2004, MNRAS, 351, 193 
Vaughan, S. & Fabian, A. C. 2004, MNRAS, 348, 1415 
van den Bergh, S., Herbst, E. & Pritchet, C. 1973, AJ, 78, 375 
Vanzella, E., et al. 2008, A&A, 478, 83 

Vasudevan, R. V., Mushotzky, R. F., Winter, L. M. & Fabian, A. 

C. 2009, MNRAS, 399, 1553 
Villforth, C, Koekemoer, A. M. & Grogin, N. A. 2010, ApJ, 723, 

737 

Walton, D. J., Roberts, T. P., Mateos, S. & Heard, V. 2011, 
MNRAS, 416, 1844 



Variability Selected LLAGN in the 4 Ms CDF-S 



Table 1 

Overview of Columns for Properties of CDF-S Galaxies 



Column Description 



1 Sequence number in the CDF-S catalog (XueetaL 2011) (i.e., XID) 

2 CDF-S name 

3 Redshift 

4 Method of redshift measurement (s = spectroscopic, p = photometric) 

5 Upper limit flag for net counts 

6 Total net counts in the 0.5-8.0 keV band or 3<r upper limit 

7 log i .5-8keV (ergs s" 1 ) 

8-10 Effective photon index (r c g ) a and corresponding errors 

11 Variability statistic (X 2 ) 

12 Probability that X 2 statistic is due to chance (P X 2 ) 

13 Maximum-to- minimum flux ratio (-R max / m i n ) 
14-15 Excess variance (cr^ xs ) b and corresponding error 

16-17 Bias-corrected excess variance (<r^ xs corr , see N4.1H and corresponding error 

18 SFR (M Q yr- 1 ) 

19 M* (10 10 M Q ) 

20 ^hmxb/^lmxb 

21 fxRB 

22 Visual galaxy classification 

23 M v 

24 Mu - My 

25 Color galaxy classification 

26 Final galaxy classification 

27 log M BH (M q ) 

28 log M/M Edd 



a The effective photon index is calculated from the band ratio. For sources detected in 
the soft band but not the hard band, T e ff is a lower limit. When the counts are too low 
to determine reliably the photon index from the band ratio, r of j is set to 1.4. 
b Erro rs on excess variance are calculated according to Equation (11) from Vaughan et al.l 
H2003T I. The excess variance is not significant for non- variable galaxies, but the errors 
can be used to calculate an upper limit. 



Table 2 

Properties of CDF-S Galaxies 



XID 


CDF-S 


name 


z a 


Net 


log 


r e ff 


X 2 




^raax / min 


2 

ff nxs 


2 

^nxs^corr 










counts 


£o.5-8koV 












120 


033206.40- 


-274728.6 


1.02 


32.7 


41.62 


1.4 


2.8 


0.0296 


4.1 


0.22±0.24 


0.56±0.63 


151 


033209.79- 


-274442.7 


0.08 


46.9 


39.39 


>0.9 


2.9 


0.0304 


4.1 


0.32±0.21 


0.85±0.75 


162 


033210.72- 


-274234.9 


0.42 


111.7 


41.16 


i o+ - 4 

1 - IJ -0.3 


5.8 


0.0002 


3.4 


0.22±0.11 


0.60±0.49 


223 


033215.80- 


-275324.7 


0.67 


<43.7 


41.51 


>1.0 


5.4 


0.0034 


8.1 


0.56±0.29 


1.67±1.34 


233 


033216.76- 


-274328.2 


0.52 


61.0 


41.26 


>1.0 


4.1 


0.0050 


3.5 


0.43±0.21 


1.16±0.97 



Note. — Table 2 is presented in its entirety in the electronic edition. A portion is shown here for guidance regarding 
its form and content. The full table contains 28 columns of information for the 20 variable and 72 non-variable CDF-S 
galaxies. 

a For 18 of 20 variable galaxies and 61 of 72 non-variable galaxies, the redshift is measured spcctroscopically and is "secure" 
(see Q. The remaining galaxies have photometric redshifts. 
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Figure 1. The fraction of variable AGN (red circles) and galaxies (black squares) vs. net counts, where the background-subtracted (net) 
counts are measured in the 0.5-8.0 keV band over 4 Ms. The binomial errors on the fraction at the 90% confidence level are calculated 
from Cameron (2011). The error bar on the net counts represents the bin size. The points are slightly offset for clarity. 
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Figure 2. The 0.5—8 keV luminosity vs. redshift distribution for CDF-S galaxies meeting the quality criteria of i|3](open circles). Galaxies 
with significant variability are marked as red stars. The CDF-S on-axis flux limit for 20 net counts (Fq.s— 8keV ~ 4.7 X 10~ 17 ergs cm -2 s _1 ) 
is plotted as a dotted line. The upper rc-axis shows the maximum rest-frame timescale in years, where the maximum observed-frame timescale 
of the CDF-S is 10.8 years. 
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Figure 3. Light curves for 6 representative CDF-S sources classified as galaxies that exhibit significant variability. Full-band (0. 5—8 keV, 
observ ed- frame) counts have been grouped into the four epochs and asymmetric errors on the count rates were calculated via IGehreld 
(1986). The mean count rate is overplotted as a dashed line, and the error on the mean is shown as the grey shaded area. Each plot lists 
the P X 2 values and source names. 
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Figure 4. The distribution of X-ray luminosity (I/o.5_ gkev) f° r an CDF-S sources classified as galaxies (dark grey histogram) and for 
those exhibiting significant variability (red histogram). The vertical lines show the median values for all galaxies (black, dotted) and for 
variable galaxies (red, dashed). A K-S test shows no significant difference between the distributions. 
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Figure 5. The fraction of variable AGN (red circles) and galaxies (black squares) vs. X-ray luminosity. Error bars show the binomial 
errors at 90% significance level. The error bar on luminosity represents the bin size. The points are slightly offset for clarity. 
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Figure 6. ij-band magnitude vs. full-band X-ray flux for sources in the main CDF-S catalog of IXue et al.| (120111 ), where filled symbols 
mark sources with significant variability (Px 2 < 0.05). Sources classified as AGN, galaxies, and stars are marked with red, black, and 
blue symbols, respectively. Constant flux ratios are marked with diagonal lines: from top to bottom, log(F x /Fn) = +1, —1, and —2. The 
vertical dotted line marks the on-axis full-band flux limit for 20 net counts. The sources with _R-band magnitude lower limits at the top of 
the plot (f»25% of CDF-S sources) have no i?-band counterpart. 
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Figure 7. The distribution of max-to-min ratios for simulated light curves based on CDF-S galaxies (black histogram) and for measured 
light curves of variable galaxies (red histogram). 
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Figure 8. The distribution of the bias-corrected normalized excess variance cr^ xs corr for CDF-S sources classified as variable galaxies (red 
histogram), with non-variable galaxies (blue histogram) shown for reference. As expected, the distribution of non-variable galaxies, though 
affected by scatter due to statistical fluctuations, centers around zero (i.e., variability strength is attributable to Poisson fluctuations). The 
red dashed line marks the maximum variability expected from an XRB population. 
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Figure 9. Postage-stamp images from the GOODS-S/GEMS HST V606 band for 20 variable galaxies. The label at the top of each 
image gives the source name. The numbers at the bottom of each image indicate the source number (XID) in the main CDF-S catalog, the 
adopted redshift, and the logarithm of the full-band luminosity as calculated in [[2] The circle overplotted on each image has a radius of 
1.5 times the X-r ay positional error, which is calculated at the 90% significance level, to illustrate whether an X-ray source is considered 
to be off-nuclear (Lchmcr ct al. 2006). Each image is 8" X 8" with the position of the X-ray source of interest located at the center. 
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Figure 10. Postage-stamp images from the GOODS-S/GEMS HST V606 band for 72 non-variable galaxies. The labels and error circles 
are the same as for Figure [9] 
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Figure 11. The bias-corrected normalized excess variance (ff„ xs corr ) vs. X-ray luminosity (Lo.5— 8kev) f° r CDF-S sources classified as 

AGN (red circles) and variable galaxies (black circles and stars) in the 4 Ms main catalog, where sources with Lo.5— 8keV>10 41 ergs s — 1 
are fitted with a weighted least-squares regression (dashed line); the shaded area shows the dispersion aroun d th e fitted line. Black stars 
highlight six variable galaxies with significant XRB contribution to the total lumi nosit y, discussed further in i|6.2l Error bars include both 
measurement errors and errors associated with red- noise scatter and sampling ( SI4.1II , Large blue squares mark the weighted means for 
each luminosity bin; error bars are the standard deviation of the data in each bin. Ultraluminous X-ray sources (ULX) are generally not 
expected above Z/q.5— 8keV ~ 10 41 ergs s _1 (vertical dotted line). 
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Figure 12. The observed 2—10 keV X-ray power (1/2— l0keV,tot) vs. that expected from an XRB population (L2— lOkeV.XRB) for variable 
(red stars) and non-variable (green squares) galaxies. Arrows designate upper limits on L2— iokcV,tot ■ CDF-S AGN (open orange circles), 
32 local galaxies (Colbert ct al. 2004), and 20 local luminous infrared galaxies ULehmer et al.ll201(ih are plo tted for reference. L ocal galaxies 
lie on ^2— lOkeV tot = ^2-lOkeV XRB ( re d solid line), with a measured dispersion of 0.34 (red dotted lines Lchmcr ct al. 2010). 
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Figure 13. The bias-corrected normalized excess variance (<r^ xs corr ) vs. accretion rate normalized by the Eddington rate (rh= M / M^dd) 
for CDF-S sources classified as AGN (red circles) and variable galaxies (black circles). A weighted least-squares regression is fitted to all 
variable sources (dashed line), and the dotted lines show the dispersion around the best-fit line. 
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Figure 14. The bias-corrected normalized excess variance (cr^ xs corr ) vs. X-ray luminosity (Lo.5— 8kov) as in Figure [TT] for sources with z 
< 1. Model relations are shown for a universal broken power-law PSD function, where the break frequency depends on Mbh and m. The 
relations are calculated for m = 9 X 10 — 3 (dashed line) and rh = 4 X 10 -4 (dash-dotted line), where Mbh = 10 5 -10 10 Mq (low to high 
£o.5-8kcv)- 



